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Abstract: Calotropis procera (Aiton) W. T. Aiton, belonging to the family Apocynaceae, is Cs evergreen plant 
species in arid and semi-arid areas of the Punjab Province, Pakistan. It grows in a variety of habitats like 
salt affected and waterlogged area, desert/semi-desert, roadside, wasteland, graveyard, forest, crop field, 
coastline, and river/canal bank. A total of 12 populations growing in different ecological regions were 
sampled to evaluate their growth, physio-biochemical, and anatomical responses to specific environmental 
condition. Population adapted to desert/semi-desert showed vigorous growth (plant height, shoot length, 
and number of leaves), enhanced photosynthetic level (chlorophyll a, chlorophyll b, carotenoids, and total 
chlorophyll), and apparent anatomical modifications such as increased stem radius, cuticle thickness, 
storage parenchyma tissues (cortex and pith), and vascular bundles in stems, while the maximum of 
midrib and lamina thickness, epidermal cells, cuticle thickness, cortical proportion, abaxial stomatal density, 
and its area in leaves. There was high plasticity in structural and functional features of these populations, 
which enable them to survive and tolerate under such hot and dry desert environment. Population of 
saline areas exhibited very critical modifications to sustain under salt prone environment. At physiological 
level, it possesses the maximum amount of organic osmolytes (glycine betaine and proline) and 
antioxidants (superoxide dismutase (SOD), catalase (CAT), and peroxidase (POD)), while at anatomical 
level, it showed intensive sclerification, large phloem region (inner and outer), pith parenchyma cells, and 
metaxylem vessels in stems and leaves. The population of dry mountains showed very distinctive features, 
such as increased shoot ionic contents (K+ and Ca?*), collenchyma and sclerenchyma thickness in stems, 
trichomes size, and numerous small stomata on abaxial surface of leaves. It is concluded that no definite 
or precise single character can be taken as a yardstick for adjudging the biomass production in this rubber 
bush weed population. 
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1 Introduction 


Calotropis procera is drought resistant and salt tolerant xerophytic C3 plant, widely growing in 
arid and semi-arid areas. It can grow in a variety of habitats like high salinity, waterlogging, 
desert/semi-desert, river/canal bank, coastline, forest, railway tract, wasteland, open and barren 
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fields, and roadside. It keeps evergreen under open field condition from young vegetative stage to 
mature reproductive stage (Boutraa, 2010; Menge et al., 2016). In the Punjab Province of Pakistan, 
the plant shows a wider geographically distribution and invades different habitats. However, the 
mechanism involving in successful distribution of this species across heterogenic environments 
are still unexplored. In this scenario, a comprehensive research was implemented to investigate 
the responses of populations specific to environmental condition. 

Rapid increase in ambient CO» concentration and temperature fluctuation of the planet is 
causing global climatic change along with significant alteration in seasonal precipitation, and 
ultimately leading to aridity towards less or semi-arid areas of the world (Koutroulis, 2019). High 
temperature and dry environment are two main factors that strongly influence the success and 
rehabilitation of plant communities. By understanding the survival mechanism of these plant 
communities, we could improve the protection of natural vegetation of these areas (Pefiuelas et al., 
2018). As in most of the terrestrial environments, water availability is considered as major 
limiting factor for photosynthesis, growth, and productivity, even in plants thrive well under arid 
conditions, and influences the population dynamics (Duan et al., 2018). 

Abiotic stresses (drought and salinity) impose harmful impacts on plant growth through various 
means such as physiological drought condition upon restricted supply of water under osmotic 
stress, disruption of cellular functions, and physiological process like photosynthesis and 
respiration through ionic imbalance or toxicity (Muhammad et al., 2021). Ultimately, it causes 
ROS (reactive oxygen species) accumulation in metabolically active sites, which may impair the 
biological structure and its process (Sachdev et al., 2021). For example, superoxide radicle, 
hydrogen peroxide, singlet oxygen, and hydroxyl radical are generated in cell, which perturb the 
biochemical events by disrupting membranous organelles and metabolites (Ahanger et al., 2017). 
In counteract of such oxidants plants produce enzymatic and non-enzymatic antioxidants. The 
most vital enzymatic antioxidants are peroxidase (POX), superoxide dismutase (SOD), and 
catalase (CAT), while non-enzymatic antioxidants included carotenoids, polyphenols, glutathione, 
ascorbate, tocopherols, anthocyanins, etc., which scavenge the most of ROS molecules (Batool et 
al., 2013). 

Seasonal falling of leaves, smallest thicker leaves, surface pubescence, organ succulence, 
epidermal and cuticle thickness, lignification of mechanical tissues, and small and fewer stomata 
that are oriented in crypts or abaxial leaf surfaces (Abd Elhalim et al., 2016), extensive root 
system, and greater water use efficiency are the main adaptations of arid plant communities 
(Fatima et al., 2021). All of these features are specific to minimize the transpirational loss, mainly 
by leaf surface, which is crucial strategy for water conservation. Reduction in vessels diameter is 
another key feature that may affect water conduction, but it minimizes the risk of embolism or 
cavitation (Kaleem and Hameed, 2021). 

C. procera (Apocynaceae) commonly known as giant milk weed, rubber bush is a perennial 
and woody, 3—5 m long evergreen shrub. It has silver green leaves, opposite and sub-sessile, 
purple-tipped flowers, and inflated pale green seed pods (Rivas et al., 2020). It is native to 
tropical and sub-tropical regions of Asia and Africa, and commonly found in the Middle East. It 
grows as wasteland plant and reproduces by seeds only. It shows vigorous growth under adverse 
climate and poor soil condition (Pompelli et al., 2019). It can withstand drought and salt affected 
environment, but the mechanism behind their success in arid and semi-arid areas is totally unclear. 
It quickly colonizes a variety of habitats like roadside, stream/river bank, lagoon edge, and 
overgrazed area via different phenomena like anemochory and zoochory. It dominates the areas of 
sandy soil with low rainfall, and considered to be indicator of degraded lands. It attracts and 
supports 80% wildlife for completion of their life cycle (Khan et al., 2019). 

Although this species has great ecological importance as it supports the several wild species 
like pollinators, i.e., honeybees, ants, butterflies, and casual visitors like wasp, moths and 
dragonfly in such open fields and deserted environments, and it might be one of major causes of 
their wide distributional ranges. But this is not enough criteria to judge the potential of species 
survival in the variety of habitats, there would be some alteration at growth, structural, and 
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functional levels in response to different environmental stresses, which determine their ecological 
success across heterogenic environments. In this scenario, this work aims to answer the following 
questions: (1) are morpho-anatomical and physiological traits of C. procera plastic in nature? (2) 
whether this species modifies its structural and functional features according to the environmental 
conditions? and (3) what are the traits that determine the ecological success of the species? This 
study was carried out to evaluate the structural and functional mechanism of different populations 
of C. procera, and relating its adaptive features to environment, which enables them to colonize a 
variety of habitats especially in arid or semi-arid areas of the Punjab Province, Pakistan. It was 
hypothesized that this species must adopt internal (anatomical) and external (morphological) 
mechanisms to cope with prevailing environmental conditions. Such adaptive components can be 
extracted and incorporated as biological markers into sensitive species to revegetate the wasteland 
of these areas. 


2 Materials and methods 


2.1 Climatic condition of the Punjab Province 


The environmental condition of the Punjab Province ranges from tropical to sub-tropical, arid to 
semi-arid, dry, and hot. Average annual precipitation of the province fluctuates from 250 to 1000 
mm, the maximum precipitation was in the monsoon season (June—August). The northern hilly 
areas receive more annual precipitation (1000 mm) as compared with southern plain areas that 
receive less than 250 mm precipitation. The maximum average temperature of a region is 48°C in 
June and the minimum is 4°C in January (Pakistan Meteorological Department). Coordinates and 
altitudes of sampling sites were recorded with the help of the global positioning system. 


2.2 Plant material selection 


Each study site was explored in an area of 1 km? for tagging the plants of each population. Six 
largest plants (replications; n=6) were subjectively selected at each collection site (habitat) 
keeping the plant distance at least 5 m. Sampling method was a single plant per replication from 
each population. Sampling points were randomly selected within each population. Sampling was 
done during August-October 2020. The youngest shoot was selected for anatomical studies. 


2.3 Study layout, soil topography, and coordinates 


This study was carried out under natural conditions (semi-arid and arid areas) in the Punjab 
Province of Pakistan, where C. procera populations were found. A total 12 populations from 
heterogenic environments were collected to evaluate their structural and functional responses (Fig. 
1). The selected regions varied by their geographical conditions, habitats, and soil types. These 
were stream/river channel (KHP: canal bank; TBJ: near Tounsa Headwork's; ADS: near Treemu 
Headwork's), mountainous range (PCK: near mountain spring; PDK: muddy dunes along roadside; 
NSR: sandy-stone hills), hyper-arid area (DIS: Thal Desert; AHP: Cholistan Desert; DGK: flat 
plain with sandy dunes) and salt affected land (LSH: saline wasteland; KBL: hyper-saline lake; 
KKL: near salt lake). 


2.4 Soil physical-chemical parameters 


The soil taken from the root zone of each population in each habitat was utilized for 
physical-chemical characteristics. Sampling was done at two depth levels, i.e., 30 and 60 cm. Soil 
pH and electrical conductivity (EC) were assessed by making saturation paste of each soil sample 
by using a combined pH and EC meter (Cond 720, WTW series; Inolab, Weilheim, Germany), 
respectively. Soil ionic contents, i.e., Na*, K^, and Ca?* were analyzed with the help of flame 
photometer (Jenway PFP-7, ELE Instrument Co. Ltd., Felsted, UK), while Cl content was 
estimated by digital chloride meter (Model 926, Sherwood Scientific Ltd., Cambridge, UK). 
Coordinates and elevation of each collection site were determined by using Google Earth Pro 
Software. 
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KHP: canal bank TBJ: near Taunsa Headwork's ADS: near Treemu PCK: near mountain 


(Rahim Yar Khan) (Muzaffargarh) Headwork's (Jhang) spring (Khushab) 
P: 102 mm; E: 76 m; P: 211 mm; E: 131 m; P: 280 mm; E: 149 m; P: 401 mm; E: 788 m; 


28?09'19"N, 70°37’51”E 30?30'46"N, 70°50’58”E 31?08'41"N, 72?08'47"E 


32?37'11"N, 72?26'50"E 


PDK: muddy dunes along NSR: sand-stone hills DIS: Thal Desert (Layyah) AHP: Cholistan Desert 


roadside (Jhelum) (Khushab) P:250 mm; E: 162m; (Bahawalpur) 
P: 152 mm; E: 204 m; P: 435 mm; E: 288 m; 31?13'41"N, 71? 48/35"E P: 141 mm; E: 136 m; 


32?35'18"N, 73*02/36"E 34?00'37"N, 71959'15"E 30?41'00"N, 71?45'33"E 


DGK: flat plain with sandy LSH: saline wasteland KBL: hyper-saline lake KKL: near salt lake 


dunes (Dera Ghazi Khan) (Bahawalpur) (Khushab) (Chakwal) 
P:236mm;E:125m; P: 140 mm; E: 127 m; P: 590 mm; E: 743 m; P: 480 mm; E: 746 m; 


30?02'56"N, 70?38'43"E 29?2756"N, 71*58/27"E 32?3718"N, 72?12'50"E 32?45'22"N, 72?42'48"E 


Fig. 1 Map showing collection sites (a) of Calotropis procera populations and its habitat characteristics (b) in 
the Punjab Province, Pakistan. P, precipitation, E, elevation. The abbreviations of these populations are the same 
in the following figures and tables. 
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2.5 Morphological parameters 


Three average sized primary shoots were clipped from the base of each plant, and average was 
used for recording morphological characteristics. Morphological parameters including plant 
height, shoot (branch) length, shoot fresh, dry weight, leaf number, and leaf area were recorded. 
Clinometer was used for plant height, while measuring tape was used for shoot length 
measurements. Shoot fresh weight (SFW) was assessed immediately after their detachment and 
then dried in oven at 65°C until constant weight to assess shoot dry weight (SDW). Finally 
reading was taken by using digital loading balance. Leaf area was measured by using the 
following formula: 

LA-LxWxCF, (1) 
where LA is the leaf area (cm?); L is the length (cm); W is the width (cm); and CF is the 
correction factor (0.75). 


2.6 Physiological parameters 
2.6.1 Ionic contents 


Plant samples were oven-dried at 65?C for analysis of inorganic ions. The ground dried samples 
were analyzed for Na*, K^, and Ca?* after H2SO4 digestion (Wolf, 1982) using a flame photometer 
(Model 410, Sherwood Scientific Ltd., Cambridge, UK). 

2.6.2 Photosynthetic pigments 

Photosynthetic pigments (chlorophyll a (chla), chlorophyll b (chld), total chlorophyll (total chl), 
and carotenoids) were assayed by using the method of Arnon (1949). Fresh leaf material (0.1 g) 
was ground in 5 mL acetone (80%) and centrifuged at 14,000 r/m for 15 min at 4°C. The 
absorbance of the filtrate was recorded at 480, 645, and 663 nm by using spectrophotometer 


(IRMECO U2020, IRMECO GmbH & Co. KG, Lütjensee, Germany). Following formulas were 
used to calculate the final value: 


Chl a =[12.7(OD663) — 2.69( OD645))x—— XW} (2) 

Chl b =[22.9(OD645) - 4 68(OD663)]X— x], (3) 

Total chl = [20.2 (OD645) - 8.02 (OD663)] x L xW], (4) 
Carotenoids = [12.7 (OD480) — 0.114(OD663)|—0.638(OD645)] / 2500, (5) 


where OD is the optical density; V is the total volume of acetone extract (mL); and W is the fresh 
weight of the sample (g). 


2.7 Biochemical parameters 
2.7.1 Total soluble proteins (TSP) 


To analyze TSP of samples, we followed the protocol of Bradford (1976). About 5 mL potassium 
phosphate buffer (pH 7.8) was used to prepare the homogenize mixture for each leaf sample. And 
0.25 g leaf sample was crushed and ground in an ice chilled pestle and mortar until the extract 
was obtained. The extract was centrifuged in ultracentrifuge machine after maintaining the 
internal temperature constant (4°C) at 12,000 r/m. Totally 5 mL Bradford reagent and 100 pL 
extract were thoroughly mixed and vortexed for 10 s. The absorbance of each sample was 
measured through UV-spectrophotometer (IRMECO U2020, Lütjensee, Germany) at 595 nm 
wavelength. 


2.7.0 Proline content 


Fresh leaf sample was homogenized in 3% sulfo-salicylic acid for estimation of free proline 
content by following Bates et al. (1973) method. The extract was subjected to incubation with 
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acid ninhydrin and glacial acetic acid at 100°C, and finally extracted by dissolving toluene. The 
absorbance of each sample containing chromophore along toluene was observed on 
UV-spectrophometer (IRMECO U2020, Lütjensee, Germany) at 520 nm wavelength. Finally, 
standard curve was established to compute the final value of free proline content on a fresh 
weight basis. 


2.7.3 Glycine betaine (GB) 


Fresh leaf sample (0.25 g) was ground in 5 mL distilled water. The extract obtained was 
centrifuged on ultracentrifugation machine at 12,000 r/m for 15 min. About 500 uL mixture 
containing 1 mL sample and 1 mL 2 N H2SOq was added in test tube and kept into refrigerator for 
90 min. Moreover, 200.0 uL potassium triiodide, 2.8 mL distilled H2O and 6.0 mL 1, 
2-dichloroethane were added (Grieve and Grattan, 1983). Finally, absorbance was recorded at 365 
nm wavelength by using UV-spectrophotometer (IRMECO U2020, Lütjensee, Germany). 


2.7.4 Enzymatic antioxidants activities 


Fresh leaf sample (0.25 g) was ground by using potassium phosphate buffer of pH 7.8 (5 mL) in 
ice chilled pestle and mortar. The extract was subjected to ultracentrifugation machine for 15 min 
at 12,000 r/m for precipitate separation. Finally, material was preserved in scientific refrigerator 
at —20?C for further analysis of enzyme activity, i.e., superoxide dismutase (SOD), catalase (CAT), 
and peroxidase (POD). 


2.71.5 SOD 


We determined SOD activity based on the principle of NBT (nitro blue tetrazolium) by following 
the method of Giannopolitis and Ries (1977). A homogenized mixture was prepared by mixing up 
of 50.0 uL sample extract, 0.1 mL L-methionine, 0.1 mL triton-X, 400.0 uL distilled water, 50.0 
uL riboflavin, 50.0 uL NBT, and 1.0 mL potassium phosphate buffer (pH 7.0) under light for 20 
min. The absorbance was recorded by using UV-spectrophotometer (IRMECO U2020, Lütjensee, 
Germany) at 560 nm wavelength. 


2.7.6 CAT 


Homogenized mixture (3.0 mL) containing potassium phosphate buffer (1.9 mL), along with 
H20» (1.0 mL), and sample extract (0.1 mL), was prepared to determine the activity of CAT 
enzyme by following Chance and Maehly (1955) method. Readings were taken using 
UV-spectrophotometer (IRMECO U2020, Lütjensee, Germany) at 240 nm wavelength. 


2.7.7 POD 


Chance and Maehly (1955) method was applied to measure POD activity. For analysis of each 
sample, 1.00 mL reaction mixture containing H202 (0.10 mL), leaf extract (0.05 mL), guaiacol 
solution (0.10 mL), and phosphate buffer (0.75 mL) was prepared. Finally, absorbance of each 
liquate was taken using UV-spectrophotometer IRMECO U2020, Lütjensee, Germany) at 470 nm 
wavelength. 


2.8 Anatomical parameters 


For anatomical studies, the collected samples from each habitat were immediately kept in 
leak-proof plastic vials containing formalin acetic-alcohol solution (v/v formalin 5%, acetic acid 
10%, ethanol 50%, and distilled water 35%), and subsequently put into acetic alcohol solution 
(v/v ethanol 75% and acetic acid 25%) after 2—3 d for longer period of preservation. A piece (1.5 
cm) of topmost internode of main branch was taken for stem anatomy and 1.5-cm piece of mature 
leaf separated from leaf-center including midrib for leaf anatomy. Free-hand sectioning technique 
was applied for permanent slides preparation. Several grades of ethanol solutions (30%, 50%, 
70%, 90%, and 100%) and xylene (100%) were used for sample dehydration and microbial 
resistance, and a combination of safranin (196) and fast green (196) for staining (Ruzin, 1999). 
Safranin stains secondary and lignified tissues (xylem vessels and sclerenchyma), while fast green 
stains primary tissues (parenchymatous tissue). After that sections were mounted under Canada 
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balsam. Micrographs of stained sections were taken on a digital camera (Nikon FDX-35, Tokyo, 
Japan) equipped with a stereo-microscope (Nikon 104, Tokyo, Japan). All parameters were 
measured through ocular micrometer calibrated with stage micrometer. 


2.9 Statistical analysis 


One-way ANOVA (analysis of variance) under CRD (complete randomized design) was applied 
on dataset for statistical analysis (Steel et al., 1997). The means were compared using LSD (least 
significance difference) by following Snedecor and Cochran (1980). PCA (principle component 
analysis) was used to check the influence of environmental factors on population's growth, 
physiology, and structural responses under various habitats. The statistical program R software 
v.4.0.1 was used to draw correlation matrix and density heatmap among various studied variables 
(species and environment) for estimation of correlation (R Development Core Team, 2017). 


3 Results 


3.4 Morphological characteristics 


All populations of C. procera responded very differently for growth attributes relevant to habitat 
types (Table 2). Plant height was higher (4.1 m) in population AHP, which was followed (3.5 m) 
by DIS population. The minimum height (1.0 m) was recorded in PDK population. Shoot length 
was the maximum in AHP population (200.4 cm), followed by DIS population (173.8 cm). It was 
the minimum (65.2 cm) in KBL population. KHP population possessed the largest leaves (213.2 
cm”), toad by that in ADS population (180.0 cm?). LSH population had the smallest leaves 
(46.8 cm‘). 

The population DIS showed the highest number of leaves (172.4), followed by the leaf number 
(155.4) recorded in ADS population. PCK population represent the smaller number of leaves 
(49.2). The maximum shoot fresh weight (664.8 g/plant) was recorded in ADS population, which 
was significantly different (633.2 g/plant) from DIS population. The minimum (410.8 g/plant) of 
this attribute was observed in KHP population. Shoot dry weight was the maximum (277.2 g/plant) 
in TBJ population, followed by that in KBL population (260.4 g/plant). It was the minimum (84.4 
g/plant) in NSR population. 


3.1.4 Shoot ionic content 


Variations regarding plant ionic contents were huge among populations (Table 2). The maximum 
value (67.6 mg/g d.wt.) of shoot Na” was observed in LSH population, followed by KBL 
population (63.4 mg/g d.wt.). Their minimum value (15.4 mg/g d.wt.) was seen in AHP and KHP 
populations. The maximum value (41.1 mg/g d.wt.) of shoot K* was noted in NSR population, 
followed by KKL population (32.7 mg/g d.wt.). Its minimum value (14.2 mg/g d.wt.) was 
recorded in DGK population. Shoot Ca?* was the maximum (28.2 mg/g d. wt.) in NSR population, 
followed by PDK population. It was the minimum (13.8 mg/g d.wt.) in DGK population. 


3.1.2 Photosynthetic pigments 


Chlorophyll a content was the maximum (1.85 mg/g f.wt.) in AHP population, followed (1.72 
mg/g f.wt.) by DGK population (Table 2). This parameter was the minimum (0.34 mg/g f.wt.) in 
LSH, KBL, and KKL populations. The maximum value (1.12 mg/g f.wt.) of chlorophyll b was 
recorded in AHP population, which was followed by that in ADS population (1.02 mg/g f.wt.). 
The minimum value (0.16 mg/g f.wt.) of this parameter was seen in three populations of LSH, 
KBL, and KKL. The population DGK showed the highest value (0.12 mg/g f.wt.) of carotenoids 
contents, this was followed (0.08 mg/g f.wt.) by AHP population. Their minimum value (0.01 
mg/g f.wt.) was noted in KKL population. Total chlorophyll content was the highest (2.97 mg/g 
f.wt.) in AHP population, followed by DGK population. The lowest value (0.50 mg/g f.wt.) was 
recorded in populations of LSH, KBL, and KKL. The chlorophyll a/b ratio was the maximum 
(3.37) in KHP population, followed by TBJ population. The minimum value (1.13) of this 
parameter was found in ADS population. The maximum value (88.0) of Chl/Caro ratio was 
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Table 2 Morphological and physio-biochemical traits of Calotropis procera populations growing in different 
habitats of the Punjab Province, Pakistan 


Stream/river channel Mountainous range Hyper-arid area Salt affected land F 
Index LSD 
KHP TBJ ADS PCK PDK NSR DIS AHP DGK LSH KBL KKL value 
Morphological trait 
at 30 25! L5 20° 10 20 3.5% 41% 24 20° L5 15° 05119'* 
pm 86.4! 944* 115.24 138.2% 150.4: 86.4. 173.8 200.4^ 100.2: 100.2; 65.2% 115.8. 153 4.8 
mes 2132* 121.64 180.00 127.2% 82.0% 112.0* 129.2% 148.8 85.2" 46.8 99.6° 146.8" 25.611.8'" 
Number of 


leaves per 108.8. 85.2 155.4 P 492? 129.60 732" 172.5 572! 62.44 102.44 85.2° 124.8° 16.110.5"" 
shoot 


Shoot fresh 
weight 410.8. 613.2 664.8* 581.65 506.49 506.49 633.2" 581.65 448.4^ 530.2! 448.4^ 550.4* 20.538.9°" 
g 
(g/plant) 
Shoot dry 
weight 181.2* 277.2 182.44 130.4! 163.2; 84.44 104.8" 130.4! 117.2% 233.2* 260.4 1632: 12.5 8.0" 
(g/plant) 
Plant ionic content 
Shoot Na* f de d d d de e f e a ab c ore 
(mggdwt) 14 29.0% 334! 322! 30.8 27.9% 2566 154" 22 67.6" 634^ 435* 7.5 20.9 
Shoot K* 20.44 19.8: 22.05 28.3% 26.9: 41.1% 23.99 19.* 142! 28.0% 264* 32.7% 53 143'" 
(mg/g d.wt.) 
Shoot Ca?* à E b " si " : b : . : : e 
(mggdwi) 199" 182" 245° 184% 267^ 282 169 242 138 168 211° 201* 34 22.2 
Photosynthetic pigment 
un T 0.91 1.22f 1.168 1.51! 1.37: 1.660 045! 1.85" 1.72 0.34 0.34 034^ 030 43.8" 
Er 0.27 0.54 1.00 0.775 0.725 0.764 0.22? 1.12? 0.865 0.165 0.16? 0.162 0.50 26.6" 
Carotenoids — 50, 0.02 0.044 0,05 0.044 0.044 0,034 0,08 0.12* 0.066 0.03% 0.01 0.02 56.27" 
(mg/g f.wt.) 
Total 


chlorophyll 1.18" 1.760 2.18 228' 2.09! 2.324 0.67 2.97 2.58" 0.50 0.50 0.50 0.40 30.377 
(mg/g £.wt.) 
Chlorophyll 
a/b ratio 


Chl/Caroratio 29.5! 88.0° 54.55 45.64 52.5" 58.0 22.38 37.5: 21.5" 83 16.6 54.0° 5.0 17.9°* 


3.87 225^" 1.3* 1.96 1.90 2.18° 2.04! 1.65f 2.019 212* 2.12* 1.705 0.20 09.5" 


Organic osmolyte 


Proline 24.99 22.5% 13.11 22.3% 13.15 19.7 22.8% 25.94 14.4 43,3" 340^ 29.5 4.5 2317" 
(umol/g f.wt.) 

Glycine 

betaine 25% 29% 14 37 34> 20 104 23° 199 49% 56 49b 15 96" 


(umol/g f.wt.) 


Total soluble 
protein (ug/g 204.08 333.7% 91.6' 208.2% 358.4° 102.1 320.8° 222.8% 247.4! 127.0" 389.15 475.6* 30.8 43.6" 


f.wt.) 


Antioxidant 
SOD . 1.02 1.4° 1.7* 135 0.9% QT 1.6° 1.1¢ 2.0° 20 lC 1.5% 03 275" 
(U/ug protein) 
CAT e e g h c f d f d b a b ney 
. 1.1 1.2 0.85 0.5 1.6 0.9'* 1.4 1.0 1.3 1.8 2.0 1.8 0.2 44.8 
(U/ug protein) 
POD 


: 1.7° 1.1° 1.43 1.0% 07 1.0% - 1.3% 177 2.18 1.5 1.6% 02 20." 
U/ug protein 

Note: Different lowercase letters within each row indicate significant difference among different populations at P«0.05 level. *, **, and 

^" indicate the significances at P«0.05, P<0.01, and P<0.001 levels, respectively. LSD, least significance difference. SOD, superoxide 

dismutase; CAT, catalase; POD; peroxidase; d.wt., dry weight; f.wt., fresh weight; - means no value. 
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recorded in TBJ population, followed by NSR population. Their minimum value (8.3) was noted 
in LSH population. 


3.1.3 Organic osmolyte 


Proline content was the highest (43.3 umol/g) in LSH population, followed by KBL population 
(Table 2). The minimum value (13.1 uimol/g) was recorded in populations of PDK and ADS. The 
maximum (5.6 umol/g f.wt.) of GB content was found in KBL population, followed by LSH and 
KKL populations. The minimum value (1.0 pmol/g f.wt) was found in DIS population. 
Population KKL showed the maximum value (475.6 ug/g f.wt.) of total soluble proteins, followed 
by KBL population. The minimum value (91.6 ug/g f.wt.) was observed in ADS population. 


3.1.4 Antioxidant 


The maximum value (2.0 U/ug protein) of SOD was recorded in populations of DGK and LSH, 
which was not significantly varied (1.7 U/pg protein) from ADS population. Their minimum 
value (0.7 U/ug protein) was noted in NSR population. KBL population possessed the highest 
value (2.0 U/ug protein) of CAT, which was close (1.8 U/ug protein) to populations KKL and 
LSH. The minimum value (0.5 U/ug protein) was found in PCK population. LSH population 
showed the maximum (2.1 U/ug protein) of POD (peroxidase), followed by KBL population. 
PDK population had a low amount (0.7 U/ug protein) of POD (Table 2). 


3.0 Anatomical characteristics 
3.2.4 Stem anatomy 


Huge variations regarding stem anatomical attributes were seen in differently adapted 
populations of C. procera (Table 3; Fig. 2). They showed remarkable modifications in size, 
nature and shape of dermal, vascular, and storage tissues. DIS population showed widened stem 
radius (1988.3 um), followed by DGK population. The narrowest stem radius (1506.7 um) was 
found in KBL population. Thick cuticle (51.5 um) was seen in AHP population, followed by 
NSR population. Thin cuticle (18.4 jum) was observed in two populations of ADS and KBL. 
NSR population showed the largest epidermal cell area (663.5 jum?), followed by DIS and KBL 
populations. PCK population indicated the smallest epidermal cell area (231.3 m7) than the rest 
of populations. 

Thick collenchyma layer (100.7 um) was recorded in PCK population, followed by DGK, LSH, 
and KKL populations. Thin collenchyma layer (48.8 um) was found in populations of AHP, NSR, 
and KHP. Enlarged cortical region (198.9 um) was noted in AHP population, followed by DGK, 
TBJ, and KHP populations. This was greatly reduced (58.9 um) in NSR population. The largest 
cortical cell area (252.2 uum?) was recorded in AHP population, followed by KHP population. 
Their smaller cell area (35.2 um?) was noted in two populations of ADS and PCK. Intensive 
sclerification (790.4 um) was seen in populations of PDK and KBL, followed by TBJ population. 
The minimum (263.3 um) was found in DGK population. 

DIS population possessed widened metaxylem area (247.0 um?), followed by LSH population. 
The narrowest xylem area (23.5 uum?) was found in KBL population. Two populations of AHP and 
LSH showed the largest outer phloem region (588.9 um’), followed by PCK population. The 
smallest region (58.3 um?) was noted in KHP population. KKL population showed the broadest 
inner phloem region (503.0 m7), its closed value was found in LSH population. The narrowest 
phloem region (85.2 um?) was noted in NSR and KHP populations. DIS population had enlarged 
pith thickness (352.1 um), followed by LSH and ADS populations. The reduced pith thickness 
(111.3 um) was found in KBL and TBJ populations. The largest pith cell area (217.2 um?) was 
found in two LSH and AHP populations, followed by KKL population. The smallest pith cell area 
(41.4 um?) was recorded in TBJ population (Table 3; Fig. 2). 
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Table 3 Stem anatomical traits of Calotropis procera populations growing in different habitats of the Punjab 
Province, Pakistan 


Stream/river channel Mountainous range Hyper-arid area Salt affected land F 


Index LSD 1 
KHP TBJ ADS PCK PDK NSR DIS AHP DGK LSH KBL KKL VUUR 


Stem radius 
(um) 
Cuticle 
thickness 364" 37.4% 18.4° 31.44 38.4% 404^ 31.4c4 51.5* 26.74 34.1° 18.4° 36.1b° 7.4 84" 
(um) 

Epidermal 

cell area 425.9% 362.7% 450.5° 231.3* 387.54 663.5 638.3* 380.09 567.3^ 348.8% 638.3% 350.2% 50.0 10.8" 
(um?) 

Collenchyma 

thickness 48.8 55.84 64.27 100.7; 55.8" 48.84 060.9" 488! 75.9 75.9% 62.6 759" 10.5 84" 
(um) 
Cortical 
region 
thickness 
(um) 
Cortical cell 
area (jum?) 


1696.7! 1633.3* 1918.3° 1696.7* 1696.7* 1955.8*^ 1988.3* 1860.0* 1918.35 1855.0* 1506.7 1849.27 59.5 2.455 


168.3^ 168.3" 83.5 $84.84 96.2° 58.94 79.3" 198.9* 168.3^ 90.0° 108.0° 106.0° 32.7 15.1" 


212.00 90.1* 35.2! 352! 105.54 153.2; 78.65 252.2? 142.1° 78.5° 99.54 143.9 33.3282" 


Sclerenchyma 
thickness 564.24 775.0° 664.6" 693.3° 790.4* 700.4^ 300.08 400.0! 263.3" 303.0* 790.4 ^ 511.0* 28.9 1.155 


(um) 


Metaxylem 


2 41.8% 105.5" 102.1] 99.0° 73.65 96.7 247.0 50.00 71.73 2012^ 23.5! 118.8" 34.128.377 
area (um^) 


Outer phloem 


2 58.3* 166.74 244.2 453.6> 105.54 136.74 229.4° 588.9* 180.6% 588.9* 118.54 214.5" 80.824.777 
area (um^) 


Inner phloem 


2 85.2" 193.24 204.24 326.2° 156.7* 85.2" 160.0% 312.9° 203.34 431.1° 116.6° 503.0* 89.120.2" 
area (um?) 


a 216.1* 111.3* 298.9% 183.34 249.8" 170.34 352.1% 217.9° 214.6° 298.9" 111.3* 186.05 49.3 14.1"* 
du als 163.2: 41.48 106.2: 148.3: 73.6" 158.5° 132.64 217.2* 12749 2172* 97.5 190.0" 31.0 1L.8" 


eee 


Note: Different lowercase letters within each row indicate significant difference among different populations at P<0.05 level. ™ and 
indicate the significances at P<0.01, and P<0.001 levels, respectively. NS, not significant; LSD, least significance difference. 


3.2.2 Leaf anatomy 


Various populations of C. procera showed very specific leaf anatomical modifications in response 
to varied environmental conditions (Table 4; Fig. 3). Two populations of DIS and KHP showed 
thick leaf midrib (5517.9 um), followed by AHP population. The minimum value (1018.3 um) 
was found in LSH population. AHP population showed very thick lamina (246.5 um), its closed 
value was seen in DGK population. Thin lamina (83.2 um) was noted in PCK population. DIS 
population possessed the largest epidermal cell area (101.0 m7), followed by LSH population. 
The smaller cell area (16.2 um?) was noted in populations of ADS and KHP. KBL and AHP 
populations had thick cuticle (35.3 um), followed by KKL and PCK populations. Thin cuticle 
(13.0 um) was noted in ADS population. Intensive sclerification (846.5 um) was observed in KBL 
population, followed by TBJ population. The reduced tissue sclerification (147.8 um) was noticed 
in ADS population. DIS population had the largest cortical thickness (1067.7 um), followed by 
KHP population. The reduced proportion (262.8 um) of cortex was found only in AHP 
population. Two populations of KKL and DIS showed the largest cells of cortex (432.0 um’), 
followed by KHP population. The smallest cells (146.6 uum?) were recorded in LSH population. 
DIS population showed widened xylem vessels (633.2 um°), followed by TBJ population. The 
narrowest xylem vessels (183.4 jum?) were found in PDK population. The maximum phloem area 
(301.2 um?) was measured in two populations of TBJ and KKL, its closed value was 


202305.00131v1 


chinaXiv 


ChinaXiv& f'ERBTIJ 


588 JOURNAL OF ARID LAND 2023 Vol. 15 No. 5 


Population Dermal tissue Vascular tissue Storage tissue 


KHP: thick cuticle (CU), 
prominent cortical parenchyma 
cells (CC), greatly enlarged 
metaxylem vessels (MV), and pith 
cell area (PC). 


TBJ: stem hairiness (SH), 
compactly arranged cortical cells 
(CC), greatly enlarged 
sclerenchyma (SC), vascular 
bundles (VB), and pith region 
(PR). 


ADS: enlarged cortical region 
(CR), metaxylem vessels (MV), 
sclerenchyma (SC), and pith 
parenchyma cells (PC) with small 
spaces. 


PCK: greater cuticle (CU), 
collenchyma thickness (CO), 
greatly enlarge cortical proportion 
(CR), metaxylem vessels (MV), 
and pith cells (PC). 


PDK: small and sparse surface 
hairiness (SH), thick cortex (CR), 
large xylem vessels (MV), 
sclerenchyma (SC), and pith 
parenchyma cells (PC). 


NSR: thick epidermis (EP), 
enlarged cortical region (CR), 
sclerenchyma (SC), inner 
phloem area (IP), and pith cells 
thin walled (PC). 


DIS: cortical cells (CC) small and 
thin walled, large metaxylem 
vessels (MV), sclerenchyma 
(SC), and loosely arranged pith 
parenchyma cells (PC). 


AHR: thick cuticle (CU), greatly 
enlarged cortical region (CR), 
metaxylem vessel (MV), inner 
phloem (IP), and pith region (PR). 


DGK: thick cuticle (CU), greatly 
reduced cortical region (CR), 
small metaxylem vessels (MV), 
sclerenchyma (SC), and pith cells 
(PC). 


LSH: stem hairiness (SH), 
compactly arranged cortical cells 
(CC), enlarged xylem vessel 
(MV), and pith parenchyma cells 
(PC). 


KBL: stem hairiness (SH), 
sclerified epidermis (EP), 
enlarged sclerenchyma bundles 
(SC), outer phloem region 

(OP), and enlarged pith cells (PC). 


KKL: thick cuticle (CU), ortical 
cells (CC) small and compactly 
arranged, enlarged sclerenchyma 
bundles (SC), inner phloem 

(IP), and pith cells (PC) small and 
loosely arranged. 


different habitats of the Punjab Province, Pakistan 


recorded in AHP population. The minimum phloem area (101.9 uum?) was observed in PCK and 
PDK populations. The largest trichome area (355.2 jum?) was found in PCK population, followed 
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by AHP population. The smallest trichome area (61.9 um?) was found only in DGK population. 
AHP population showed very dense trichome on leaf surface (119.2/mm?), followed by PCK 
population. The lowest density (40.5/mm?) was noticed in two populations of DGK and KBL 
(Table 4; Fig. 3). 


Table 4 Leaf anatomical traits of Calotropis procera populations growing in different habitats of the Punjab 
Province, Pakistan 
Stream/river channel Mountainous range Hyper-arid area Salt affected land 


Index LSD F value 
KHP TBJ ADS PCK PDK NSR DIS AHP DGK [LSH KBL KKL 


Midrib 
thickness 5517.9? 4273.7° 2288.32 3365.6! 4132.14 4135.44 5517.97 4398.3^ 1198.3" 1018.3! 3796.2* 1195.0^ 500.9 272.3" 
(um) 

Lamina 

thickness 176.1° 103.9% 198.9° 83.25 223.6" 176.4° 183.34 246.5* 232.1% 194.7 132.0! 133.8' 11.2 80.4" 
(um) 

Epidermal 

cell area 16.22 39.3° 16.2% 20.1! 85.0" 21.6! 101.0* 85.9b° 74.3: 89.35 46.00 84.35 4.7 201.5"** 
(um?) 

Cuticle 

thickness 28.8 32.0 13.05 33.4% 19.34 16.0 18.75 35.3 28.0 18.84 35.3 33.4% 2.7 343'" 
(um) 

Scleren- 

chyma 723.0° 771.3 147.8) 155.4! 296.8" 504.3* 501.5° 154.5' 448.1 313.0? 846.5* 548.44 39.4 104.27" 
thickness 

(um) 

Cortical 

thickness 990.5> 641.67 896.6° 389.9% 614.0% 610.4% 1067.7? 262.8" 672.44 853.0° 567.8° 470.4' 75.8 42.3" 
(um) 

Cortical 

cell area 385.3? 251.2: 161.3" 231.4! 176.3% 260.7° 432.07 294.34 176.5? 146.6! 307.7° 432.0? 21.3 530.6°™ 
(um?) 
Metaxylem 
area (um?) 
Phloem 
area (um?) 


530.1^ 608.5? 511.5" 423.4° 183.4! 266.1* 633.2* 583.69 443.2 364.34 537.7” 493.8" 76.1 19.0" 


139.45 301.2* 161.3! 101.9! 101.9! 232.94 166.7! 285.3 115.6" 172.6* 242.7 301.2? 232 303.3" 


Trichome 


>, 137.8" 261.4 207.3* 355.27 108.1! 144.4? 81.2) 280.1” 61.9% 258.14 174.9 81.7) 19.8 929.3" 
area (um*) 


Trichome 
density 


f d b e a h c h c Ett 
(individu- 52.6, 79.9" 49.38 100.8^ 49.98  69.3* 49.5% 119.2? 40.5 85.3° 40.5 90.3* 6.6 151.0 


als/mm?) 


Adaxial 

stomatal 914.2 827.7? 957.3* 432.2) 700.8% 685.1^ 682.6" 796.0* 718.0! 496.3! 795.5° 861.4° 21.0 1657.5" 
area (jum?) 

Abaxial 

stomatal 842.9? 524.45 737.00 100.3! 83.8 100.2! 269.0% 881.4, 83.8) 121.5" 508.6° 404.9! 20.5 3545.1" 
area (um?) 

Adaxial 

stomatal 

density 35.00 50.0 * 45.0" 40.0" 35.0? 32.0% 32.0% 23.05 29.0 " 26.00 231.0 28.0 5.8 100.2" 
(individu- 

als/mm?) 

Abaxial 

stomatal 

density 27.0° 26.00 29.0" 21.00 31.0 16.00 20.0 31.0 27.0 21.0° 29.0% 29.0% 232 73.9 
(individu- 

als/mm?) 


ok 


Note: Different lowercase letters within each row indicate significant difference among different populations at P«0.05 level. 
indicate the significance at P<0.001 level. LSD, least significance difference. 
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Population Dermal tissue Vascular tissue Storage tissue 
Pa M = 


KHP: greatly reduced epidermal 
thickness (EP), metaxylem 
vessel area (MV), enlarged outer 
phloem (OP), and cortical region 
(CR). 


TBJ: thick cuticle (CU), enlarged 
collenchyma bundles (CO), xylem 
vessels (MV), outer phloem area 
(OP), and cortical parenchyma 
cells (CC). 


ADS: thick cuticle (CU), 
collenchyma (CO), xylem vessels 
(MV) few and small, and enlarged 
cortical parenchyma cells (CC). 


PCK: thick cuticle (CU), 
increased collenchyma thickness 
(CO), greatly reduced vascular 
bundles (VB), and cortical cells 
large (CC). 


PDK: thick epidermis (EP), 
reduced metaxylem vessels (MV), 
and irregularly arranged cortical 
parenchyma cells (CC). 


NSR: enhanced epidermal 
thickness (EP), collenchyma layer 
(CO), enlarged vascular bundles 
(VB), and cortical parenchyma 
cells (CC). 


DIS: greatly enlarged epidermal 
cells (EP), xylem vessels (MV), 
outer phloem region (OP), and 

cortical parenchyma cells (CC). 


AHP: thick cuticle (CU), well 
developed collenchyma (CO), 
metaxylem vessels (MV), 

and small and loosely arranged 
cortical cells (CC). 


DGK: increased cuticle thickness 
(CU), collenchyma thickness 
(CO), metaxylem area (MV), 
outer phloem region (OP), and 
enlarged cortical region (CR). 


LSH: reduced epidermal thickness 
(EP), vascular region (VB), and 
cortical parenchyma disintegrated 
(CR). 


KBL: increased cuticle thickness 
(CU), collenchyma bundles (CO), 
greatly reduced vascular system 
(VB), and enlarged cortical 
parenchyma cells (CC). 


KKL: thick cuticle (CU), enlarged 
collenchyma region (CO), 
metaxylem vessels (MV), and 
cortical cells (CC) irregularly 
arranged. 


Fig. 3 Transverse sections showing different leaf tissues of Calotropis procera populations collected from 
different habitats of the Punjab Province, Pakistan 


The maximum stomatal size (957.3 um?) on adaxial surface of leaf was noticed in ADS 
population, followed by KHP population. The minimum size of stomata (432.2 um?) was 
recorded in PCK population. On the other hand, the largest stomatal size (881.4 uum?) on abaxial 
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surface was measured in AHP population, followed by KHP population. The smallest size of 
stomata (83.8 um?) was recorded in two populations of PDK and PCK. The adaxial stomatal 
density was the maximum (50.0/mm?) in TBJ population, followed by ADS population. The 
minimum stomatal density (23.0/mm?) was found only in AHP population. The maximum abaxial 
stomatal density (31.0/mm?) was recorded in two populations of AHP and PDK. Its closed value 
was seemed in three populations of KBL, KKL, and ADS. The minimum of this attribute 
(16.0/mm?) was observed in NSR population (Table 4; Fig. 4). 


3.3 Multivariate analysis 


PCA was applied on data to check the influence of environmental factors on growth, 
physiological, and anatomical attributes of different populations of C. procera (Fig. 5a-d). 

3.3.1 Influence of soil physical-chemical attributes on growth and physiology 

Plant height, leaf number, and shoot length were strongly influenced in AHP population, while 
leaf area was greatly influenced in KHP population rich in soil Ca?*. Shoot dry weight was 
influenced by soil EC, Na‘, and Cl’ of LSH and KBL populations, whereas shoot fresh weight 
was greatly influenced by soil K* in DIS population. In case of plant physiology, SOD, shoot K^, 
shoot Na*, and proline showed strong impacts of soil Na* in KKL population, while GB, CAT, 
and POD represent the effects of soil EC and Cl in KBL population. Photosynthetic pigments 
like Chlorophyll a, Chlorophyll 5, and total chlorophyll showed strong association with AHP 
population rich in soil K^, but the carotenoids showed association with PCK population 
impregnated by soil Ca?*. 

3.3.2 Influence of soil physical-chemical attributes on anatomy 


In case of stem attributes, sclerenchyma thickness showed positive relation with KBL population, 
while collenchyma thickness showed positive relation with LSH population under varying level 
of soil EC, Na, and CI. Cuticle thickness and stem radius showed association with DGK 
population, whereas cortical cell area, cortical region thickness, and epidermal cell area showed 
associated with soil pH in DIS population. Among leaf anatomical attributes, cuticle thickness 
was strongly influenced in KKL population, while trichome area and trichome density were 
influenced by soil Ca?* in PCK population. In NSR population, metaxylem area and lamina 
thickness were associated with soil K^, adaxial stomatal area and abaxial stomatal area in TBJ 
population, while adaxial stomatal density and cortical region thickness showed close association 
with AHP population. 


3.4 Correlation matrix and density heatmaps 


The soil physical-chemical traits and morpho-physiological traits exhibited a significant (P«0.01) 
correlation (Fig. 6a). Soil physical-chemical traits, i.e., Cl, Na‘, and EC reflected a strong 
positive correlation with osmolytes like proline and GB, and with antioxidant enzymes (CAT and 
POD). Similarly shoot Na was positively correlated with organic osmolytes and cellular 
antioxidants. Soil and shoot K* and Ca?* were negatively correlated with antioxidants. Clustered 
heatmap was constructed between soil physical-chemical traits and morpho-physiological traits to 
visualize the association between studied traits. Soil EC, Na*, and Cl” exhibited a strong relation 
with shoot Na‘, shoot dry weight, POD, and SOD, while a negative association was assessed 
between photosynthetic pigments and soil traits. Moreover, a distinct clustering was observed 
between soil traits and morpho-physiological characters (Fig. 6b). 

The soil physical-chemical characters possessed a significant influence on stem structural 
attributes as presented in Figure 7a. The positive correlation was passed between Na” and inner 
phloem area, stem radius, pith thickness, inner phloem area, and pith cell area. However, a strong 
negative influence was observed of Ca?' with epidermal cell area, pith cell area and 
sclerenchymatous thickness. Soil Cl, Na*, and EC indicate a strong association with outer 
phloem area, inner phloem area, and cortical cell area. A negative association was observed 
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Population 


KHP: large stomata (ST), few 
with 4—5 irregular size 


Ü 
subsidiary cells (SC), and small F 


surface hairiness (SH). 


TBJ: large stomata (ST), 
numerous with irregular 


subsidiary cells (SC), small and = z : ^ 


sparsely arranged surface 
hairiness (SH). 


ADS: large stomata (ST), 
bulbous and numerous with 4— 
5 subsidiary cells (SC), and 
small densely arranged surface 
hairiness (SH). 


PCK: small stomata (ST), 
numerous and elliptic with 
irregular subsidiary cells (SC), 
and few large trichomes (TR) 
are sparsely arranged. 


PDK: large stomata (ST), 
numerous with 4—5 

subsidiary cells (SC), and 
surface hairiness (SH) is much 
reduced. 


NSR: small stomata (ST), 
sunken and numerous with few 


subsidiary cells (SC), and small @& 1 


surface hairiness (SH). 


DIS: large stomata (ST), few 
with reduced subsidiary cells 
(SC), and surface hairiness 
(SH) is reduced and 

sparsely arranged. 


AHP: large stomata (ST), few 
on adaxial side, while small on 
abaxial side, and surface 
hairiness (SH) is very sparse. 


DGK: large stomata (ST), few 
surrounded with 4—5 subsidiary 


cells (SC), and surface hairiness f 


(SH) is very reduced. 


LSH: large stomata (ST), 
numerous with reduced 
subsidiary cells (SC), 

and surface hairiness (SH) is 
very small. 


KBL: small stomata (ST), few 
arranged with irregularly 
subsidiary cells (SC), and surfac 
hairiness (SH) is very sparse . 


KKL: numerous small stomata 
(ST), with more than 4 
subsidiary cells (SC), and few 
trichomes (TR) are 

sparsely arranged. 


Abaxial stomata 


Surface view 


Fig. 4 Transverse sections showing epidermal surface view of Calotropis procera populations collected from 


different habitats of the Punjab Province, Pakistan 
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Fig. 5 PCA (principle component analysis) showing association of soil physical-chemical characteristics with 
growth (a), physiology (b), stem (c), and leaf anatomical attributes (d) of Calotropis procera populations collected 
from different habitats of the Punjab Province, Pakistan. PH, plant height; LN, number of leaves; LA, leaf area; 
SL, shoot length; SFW, shoot fresh weight; SDW, shoot dry weight; SNa*, shoot Na*; SK*, shoot K*; SCa?*, shoot 
Ca?*; Chla, Chlorophyll a; Chlb, Chlorophyll b; Caro, carotenoids; Tchl, total chlorophyll; Pro, proline; GB, 
glycine betaine; TSP, total soluble proteins; SOD, superoxide dismutase; CAT, catalase; POD, peroxidase; Ca?', 
soil Ca?*; K*, soil K*; Na‘, soil Nat; CL, soil CI; EC, electrochemical conductivity; pH, soil pH; CuT, cuticle 
thickness; StR, stem radius; CCA, cortical cell area; EpCA, epidermal cell area; CRT, cortical thickness; SCT, 
sclerenchymatous thickness; ColT, collenchymatous thickness; IPA, inner phloem area; OPA, outer phloem area; 
PCA', epidermal cell area; MtA, metaxylem area; PhlA, phloem area; TA, trichome area; TD, trichome density; 
LmT, lamina thickness; AdSD, adaxial stomatal density; CorT, cortical region thickness; AdSA, adaxial stomatal 
area; AbSA, abaxial stomatal area; AbSD, abaxial stomatal density; MrbT, midrib thickness. The abbreviations 
are the same in figures 6, 7 and 8. 


between epidermal cell area, cortical cell area, outer phloem area, inner phloem area, and stem 
radius with pith thickness (Fig. 7b). A correlation matrix between soil physical-chemical 
attributes and leaf structural traits indicated a significant (P<0.01) correlation. A strong positive 
correlation was assessed of trichome area and trichome density with Ca?*. The phloem area, Ca?*, 
Na‘, and adaxial stomatal density presented a significant negative correlation (Fig. 8a). Leaf 
anatomical traits exhibited a strong clustering with soil K* and Ca‘, while sclerenchymatous 
thickness grouped with midrib thickness, epidermal cell area with lamina thickness, and abaxial 
stomatal area with adaxial stomatal area. Soil Cl, Na‘, and EC showed a strong influence on 
phloem area, cortical cell area, abaxial stomatal area, adaxial stomatal area, and abaxial stomatal 
density. Density pattern was varied among all studied traits (Fig. 8b). 
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Fig. 6 (a), correlation among soil physical-chemical traits and morpho-physiological attributes, and (b), density 
heatmap for soil physical-chemical traits and morpho-physiological attributes of Calotropis procera populations 
collected from different habitats of the Punjab Province, Pakistan. “, P«0.05 level; ***, P«0.01 level. 


4 Discussion 


Huge variations regarding growth attributes were recorded among populations due to great 
environmental heterogeneity. These characteristics were genetically fixed over prolonged time of 
evolution in these populations, and these are representative of their respective habitats (Paccard et 
al., 2013). Growth attributes are considered as important criteria to estimate the response and 
tolerance potential of population to prevailing stressful condition (El-Hendawy et al., 2017). In 
this study, the individuals adapted to hyper-arid areas like AHP and DIS had the highest height 
and shoot size, compared with individuals of other habitats. It may be apparently due to 
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Fig. 7 (a), correlation among soil physical-chemical traits and stem structural attributes; (b), cluster heatmap for 
soil physical-chemical traits and stem structural attributes of Calotropis procera populations collected from 
different habitats of the Punjab Province, Pakistan. *, P<0.05 level; **, P<0.01 level; ***, P«0.001 level. 


less anthropogenic activities (i.e., cultivation, cutting, and burning), reduced grazing impact, and 
episodic water resources like annual precipitation and flooded irrigation. In earlier reports, it was 
clearly mentioned that Calotropis is a drought resistance and salt-tolerant plant, and can grow in 
different ranges of soil type. The populations growing in salt affected lands showed much reduced 
growth attributes such as plant height, shoot growth, leaf number and size, and biomass 
production. It could be result of depressive effect of salinity on growth attributes. Many researchers 
have reported that it can tolerate salinity to some extent and cannot sustain under high soil salinity 
(Boutraa, 2010; Hassan et al., 2015). 

Leaves are main photosynthetic organs of plants and thus exceptionally vital for growth and 
biomass accumulation (De Faria et al., 2018). Formation of large leaves in KHP is a principal 
strategy in order to prevent from seasonally occurring water shortage due to low precipitation 
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Fig. 8 (a), correlation among soil physical-chemical traits and leaf structural attributes; (b), cluster heatmap for 
soil physical-chemical traits and leaf structural attributes of Calotropis procera populations collected from 
different habitats of the Punjab Province, Pakistan. *, P<0.05 level; ", P«0.01 level; ***, P«0.001 level. 


and canal irrigation. According to Tezara et al. (2011) and Rivas et al. (2020), plants showed high 
photosynthetic rate, stomatal control, and efficient metabolic performance under drought stress 
condition. Leaves per plant were exceptionally high in DIS, but smaller in size. Large number of 
narrow leaves may aid in carbon acquisition as well as water use efficiency by controlling 
transpiration rate (Mutwakil et al., 2017). Fresh and dry biomass production is a good criterion to 
judge the tolerance potential in plants against abiotic stresses (Zokaee-Khosroshahi et al., 2014). 
ADS showed the maximum shoot fresh weight, while that TBJ showed the maximum of shoot dry 
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weight. The robust growth of both populations is a reflection of their complete dominance to 
respective habitats. 

Abiotic stresses like salinity and water logging, high temperature, and irradiance impose 
substantial changes in tissues ionic contents, photosynthetic pigments, organic osmolyte level, 
and antioxidants in plants (Hasanuzzaman et al., 2020). Accumulation of high Na* concentration 
in tissues has been reported in various glycophytes and halophytes (Sekmen et al., 2013), as 
recorded in LSH population in this study. These plants utilize such ions to maintain the cell 
osmotic potential and osmoregulatory process under osmotic stress condition (Rouached et al., 
2010). Remarkable concentration of shoot K^ and Ca? in NSR population showed soil 
enrichment of these ions. In addition, these are beneficial for plants to regulate the cellular 
osmotica and metabolic process (Alzahrani et al., 2019). Much higher concentration of 
photosynthetic pigments (chla, chlb, and total chl) was found in AHP population, whereas large 
quantity of carotenoids was found in DGK population. This increase in pigments seems to be 
directly associated with photosynthetic efficiency of species under stressful environment 
(Maimaitiyiming et al., 2017). 

Salt tolerant species accumulate bulk amount of organic osmolyte to offset the effects of ROS 
and free radicals on photosynthetic apparatus and organelles involved in respiratory process, that 
ultimately lead to the oxidative stress. Osmoprotectants also maintain the protein structure, 
membrane integrity, and osmotic events in plants (Muchate et al., 2016). In our results, the 
populations from salt affected lands (LSH, KBL, and KKL) were rated as more tolerant due to the 
maximum accumulation of proline, GB, and soluble proteins. Plants generally evolve antioxidants 
mechanism (enzymatic and non-enzymatic) to mitigate and detoxify the effects of ROS. 
Enzymatic antioxidants like SOD, CAT, and POD, are designed in chain reaction to eliminate and 
scavenge the superoxide and H2O» in plants (Sheikh-Mohamadi et al., 2017). In current study, all 
the populations showed different level of antioxidant metabolism with respect to prevailing stress 
conditions, while the population inhabiting in saline environment (LSH, KBL, and KKL) 
surpassed in antioxidant activity. Keeping the high level of antioxidants might supplement for 
salinity or drought stress tolerance by improving the protective strategy against oxidative stress 
(Zhang et al., 2013). 

Hameed et al. (2012) found that desert plants were adapted to limited water availability, 
therefore, the impact of water stress was not severe as on the other habitat with high moisture 
availability, especially in tolerant species. The response of differently adapted populations of C. 
procera was very specific, indicating their capability to survive in desert environment. All the 
populations showed great diversity in stem anatomical tissues like dermal, mechanical, vascular, 
and storage tissues. The population growing in desert and semi-desert environment showed the 
maximum stem anatomical parameters, i.e., cuticle thickness, cortical region thickness, cortical 
cell area, outer phloem area, and pith cell area in AHP population; and stem radius, pith region 
thickness, and metaxylem area in DIS population. Desert plants modify their structural and 
functional features to cope with multiple abiotic stresses and successfully survive in such 
voracious climatic conditions. In earlier studies, cuticle or epidermis thickness is reported as first 
line of defense in arid plants to minimize the water loss under water scarce conditions (Liu et al., 
2020). Another adaptation in these plants is the formation of enlarged stem area due to deposition 
of storage parenchyma and large vascular bundles. This may certainly aid in water conservation in 
plant tissues and hence survive with subject to prolonged dry weather condition and low moisture 
contents (Su et al., 2019). Therefore, AHP population showed enormous succulence in form of 
storage parenchyma (cortex) to survive under such environment. The formation of large phloem 
region is another principle strategy of this population that may enhance photosynthetic 
translocation and partitioning to other parts (Hussain et al., 2021). 

Vascular bundles with large metaxylem vessels play a key role in efficient translocation of 
water and minerals. Tissue sclerification is a prominent feature of plants surviving under saline 
and drought stress conditions (Hameed et al., 2020). Same result was observed in stem of KBL 
and NSR populations. Tissue sclerification gives mechanical strength to metabolically active 
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tissues in order to protect from collapse and desiccations, when water has become vital 
commodity for plants either due to low soil water potential or low water table (Abd Elhalim et al., 
2016). Thick collenchyma layer formed under epidermis and around the cortical tissues is very 
critical feature of PCK population, which might be their habitat friendly behavior or strategy to 
overcome the prevailing stress. Collenchyma acts as supporting tissue in plants of dry 
environment by preventing water loss from surface (Iqbal et al., 2021). 

Leaves are quite responsive to climatic conditions, because they are directly exposed to 
environmental stress (Drake et al., 2019). Leaf thickness showed very differential response to 
varied environmental conditions, i.e., the populations growing in desert regions showed very thick 
leaves regarding lamina (AHP) and midrib thickness (DIS and KHP). Leaf thickness might be the 
result of increased storage parenchyma region to store additional water. This is a critical adaption 
of desert dwelling species to cope with severe drought condition (Grubb et al., 2015). In addition 
to thick epidermis and cuticle, trichome size and number, stomatal size, shape, orientation, and 
regulation are also vital for plants to survive under dry environment. Stomata are generally 
oriented on abaxial side of leaves in most of desert species, to reduce transpiration rate, especially 
when leaves are directly subjected to sun rays (Silva et al., 2014). Similar result was observed in 
AHP population, having more stomata on leaf abaxial surface of plants. Moreover, small size 
stomata show efficient translocation as compared with larger ones due to low need of water for 
turgor maintenance (Rudall et al., 2013). 

Intensive sclerification around cortical and vascular region of KBL population, hence making 
leaves stiffer and more fibrous. Cuticle thickening in this population additionally contributes in 
leaf stiffness, moreover, it also protects from solar irradiance and pathogens infection (Akcin et 
al., 2015). Large metaxylem formed in leaves of KKL population might be adaptability in order to 
maintain water supply under low soil water potential. Trichomes act as insulating agent to high 
wind pressure, rain splashes, and high temperature. These are ecologically important for species 
of dry environment, which maintain the maximum humidity on surface of leaves (Dolatabadian et 
al., 2011). Few and large trichomes in PCK population indicate its tolerance and strategy to make 
sure its survival under limited water supply. ADS population growing along the river/canal bank 
possessed large stomata on adaxial leaf surface. As stomata are vital for transpiration and gaseous 
exchange in plants, its opening allow the maximum water utilization along with carbon 
sequestration. They maintain the plant's vigor and vitality by producing more energy in saturated 
environmental conditions (Buckley, 2019). 


5 Conclusions 


It is concluded that survival responses in C. procera populations appeared as the result of 
complex and differential reactions to prevailing stress conditions. The maximum modulation of 
growth, photosynthetic, and anatomical features in DIS and AHP populations strongly favored in 
their life cycle completion, and also indicated its tolerance potential under harsh environment. 
The population inhabiting saline areas showed enhanced level of metabolites and antioxidant to 
cope with osmotic stress condition. Accumulation of beneficial ions (K^ and Ca?) in NSR 
population and the maximum biomass production (shoot fresh and dry weight) in ADS population 
shows their success in respective habitats. Overall, growth, physio-biochemical, and anatomical 
indices observed in these populations are best indicators of their survival and stress tolerance. 
However, there is a need for further study into the effects of habitats on functional attributes of C. 
procera populations at molecular and genetic levels. 
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